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Abstract. Spectral datasets obtained at radio frequencies and optical/IR 
wavelengths are increasing in complexity as new facilities and instruments come 
online, resulting in an increased need to visualize and quantitatively analyze 
the velocity structures. As the visible structure in spectral data cubes is not 
purely spatial, additional insight is required to relate structures in 2D space 
plus line-of-sight velocity to their true three-dimensional (3D) structures. This 
can be achieved through the use of models that are converted to velocity-space 
representations. We have used the S2PLOT programming library to enable 
intuitive, interactive comparison between 3D models and spectral data, with 
potential for improved understanding of the spatial configurations. We also 
report on the use of 3D Cartesian shapelets to support quantitative analysis. 



1. Interactive Visualization of Spectral Data Cubes 

Spectral data cubes obtained at radio frequencies (e.g. via aperture synthesis, 
mosiacs or multibeam instruments) and optical/infrared wavelengths (e.g. via 
integral field units, scanning Fabrey-Perot interferometers, or multi-object spec- 
trographs) comprise two spatial dimensions and one spectral dimension. In both 
regimes, the spectral dimension is a proxy for the relative line-of-sight velocity 
between the observer and source, resulting in spectral features that are Doppler- 
shifted from their rest frequency. The challenge is to intuitively understand the 
relationship between the two spatial plus one velocity observations, and the true 
three-dimensional spatial plus three-dimensional velocity structure. 

Solutions exist for simple velocity structures, such as rotating extragalactic 
HI disks, which can be fit by a set of inclined, differentially rotating annuli as 
a function of radius (e.g. Rogstad et al. 1974, Begeman 1989). These models, 
however, do not easily account for kinematic features such as warps, anomalous 
gas, or mergers, and are not appropriate for interpreting the complex velocity 
structures that occur in observations of galactic-plane neutral hydrogen. 



1 



2 



Fluke, English, Barnes 




Figure 1. Mapping from simulation space (left) to spectral cube space 
(right) using the orientation of a virtual camera to determine the mapping. 
Data is from an TV-body simulation of a galaxy merger. The spectral cube 
has been oriented independently to show velocity-space features. 



One option that presents itself is to generate particle models, where the 
full 3D spatial, x = (x,y,z), and 3D velocity, v = (v x ,v y ,v z ), information is 
known for each particle, and project this to spectral cube "space", (x',y',v' z ), 
for display. The orientation of the model data with respect to the world camera, 
with unit up vector, u = (u x ,u y , u z ), view direction, w = (w x ,w y , w z ), and right 
vector, r = w x u = (r x ,r y ,r z ), is used to obtain the coordinates for plotting 
points in spectral cube space: 

x' = xr x +yr y + zr z , (1) 
y' = xu x + yu y + zu z , and (2) 

V 'z = -{V X WX + VyWy + V Z W Z ) + V sys , (3) 

where v sys is the overall system velocity. It is straightforward to convert v' z to a 
frequency or wavelength. The result of this process is demonstrated in Figure Q] 
for a galaxy merger. 

We are using the S2PLOT 3D-graphics library (Barnes et al. 2006) to de- 
velop a framework where models can be examined interactively alongside their 
spectral-cube projections. This approach enables an increased intuitive under- 
standing of complex kinematical structures. As a simple, yet powerful appli- 
cation development library (C/C+- (-/Fortran), S2PLOT provides customizable 
user interaction via keyboard and mouse controls; output to mono- and stereo- 
scopic displays; and support for publication in 3D-PDF format (e.g. Barnes 
&; Fluke 2008) via VRML-export with a single key press. Simulation inputs 
and spectral cube data formats can be handled with custom code or standard 
C libraries (e.g. cfitiso for FITS files). Figure [2] shows a screenshot from a 
prototype application - screen elements are described in the caption, and image 
annotations identify the main S2PLOT functionality we use. 

We note our S2PLOT approach is limited by the workstation's graphics 
memory, and may not be suitable for very large spectral data cubes (> 512 3 
voxels). To this end, we anticipate integrating our approach with a graphics 
processing unit (GPU) based volume rendering system (see Hassan et al. 2009). 
Alternatively, it may be sufficient to select sub-volumes to analyze. 
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Figure 2. Screenshot from a prototype S2PLOT application for interactive 
visualization of model and spectral cube data. The top left panel is model 
space, the top right panel is spectral cube space, and the bottom panel is an 
interactive rotation curve. Annotation on the figure gives the main S2PLOT 
functionality we are using. The model in this example is a rotating disk. 

We intend to include support for the following basic functionality: gener- 
ation of standard model components such as rotating disks, warps, hot-spots, 
expanding shells, iV-body simulated data; model-specific controls (e.g. interac- 
tive rotation curves and density profiles for disk models); realistic noise; overlay 
of simulated spectral data on observational data; and integrated quantitative 
analysis tools, such as 3D Cartesian shapelet decomposition. 

2. Quantitative Analysis with 3D Cartesian Shapelets 

The one- and two-dimensional Cartesian shapelets formalism was introduced by 
Refregier (2003). In one dimension, shapelets form an orthonormal basis set: 



where n is the shapelet order, H n (x), is the n-th Hermite polynomial, and (3 is a 
scale factor. We have recently extended the Cartesian approach to 3D, deriving 
a number of important analytic relationships (Fluke et al. in prep). While a full 
presentation is beyond the scope of this paper, the orthonormality of shapelet 
states enables us to write 3D shapelets in terms of functions: 



where n = (ni, 712,^.3) are the integer shapelet orders. An arbitrary (sufficiently 
well-behaved) 3D structure, /s(x), can be decomposed into a weighted sum of 



4> n {x) = {2 n ^/ 2 n)r 1/2 H n {x)e- x2 l 2 
B n (x;P) = /r 1/2 <M/3-^) 



(4) 
(5) 



(6) 
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Figure 3. Shapelet decomposition of a spectral data cube with 57 3 voxels. 
(Left) input cube; (middle) shapelet reconstruction; (right) residual cube. 
The basic features of the input cube are well-recovered. 



shapelet coeffecients, ,f 3in , by integration over a volume, V: 

f3,n= [ / 3 (x)£ 3 ,n(x;/?)d 3 X (7) 

Jv 

subject to m + ri2 + n 3 < n max , and an "optimal" /3-value. The shapelet recon- 
struction is: 

/ 3 (x)= Yl «(n)/ 3 , n S 3in (x;/3), (8) 

ni,n2,n 3 

where a(n) is an optional filter term. An example of the decomposition and re- 
construction process for a mock spectral data cube (no noise) is shown in Figure 
[3l As a highly parallel task, the decomposition and reconstruction processes can 
be efficiently implemented on a GPU (Barsdell et al. 2009). 

3. Conclusion 

We are using the S2PLOT programming library to enable interactive, visualization- 
directed model-fitting of spectral data cubes, leading to improved intuitive un- 
derstanding of complex kinematic structures. Shapelet space provides new op- 
portunities for quantitative analysis of complex 3D structures. We intend a 
public code release on completion. 
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